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c/3 , A mechanism of stimulated emission of electromagnetic radiation by an electron beam in carbon 

nanotubes is theoretically considered. Three basic properties of carbon nanotubes, a strong slowing 

C^ ■ down of surface electromagnetic waves, ballisticity of the electron motion over typical nanotube 

length, and extremely high electron current density reachable in nanotubes, allow proposing them 
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^ ' wave tube and free electron laser. Dispersion equations of the electron beam instability and the 

threshold conditions of the stimulated emission have been derived and analyzed, demonstrating 
realizability of the nanotube-based nanoFEL at realistic parameters of nanotubes and electronic 
^^ ' beams. 
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I. INTRODUCTION 

Since the discovery of carbon nanotubes (CNTs) in 1991-, there has been great interest to 
their outstanding structural, electrical and mechanical properties^i^ due to wide applications 
ranging from chemical and biological sensors and actuators to field emitters to mechanical 
fillers for composite materials. Among others, the study of CNTs as building blocks for 
nanoelectronics^ and nanooptics^ has continued to grow unabated owing to the great po- 
tentiality for the miniaturization and the increase of operational speed of optoelectronic 
nanocircuits, and for the use in near-field subwavelength optical element. In that relation, 
the question of electromagnetic response properties of CNTs arises. Many interesting phys- 
ical effects have been revealed, such as excitation of surface plasmons^, guiding of strongly 
slowed-down electromagnetic surface waves^*^, antenna effect - controlled and enhanced radi- 
ation efficiency in infrared and terahertz range^^^^^^^^^^^^, enhanced spontaneous decay rate 
of an excited atom in the vicinity of CNT— , and formation of the discrete spectrum in ther- 
mal radiation of finite-length metallic CNTs in the terahertz range^-. Recently, nanoscale 
optical imaging of single-walled CNTs has been studied by means of high-resolution near- 
field Raman microscopy^ii^ and antenna operation of a CNT array has been demonstrated 



experimentally^^. Reference |l9| reports multi-wall CNT as subwavelength coaxial waveguide 
for visible light. 

An intriguing problem of nanoelectromagnetism is the development of CNT-based 
nanoscale sources of light. A mechanism of the emission of hard X-radiation by a charged 
particle moving in a CNT has been considered in Ref. |20|. The use of CNTs in X-ray 
and high-energy particles optics as focusing and guiding elements and as x-rays sources is 
presently discussed^i^. In the optical range, the mechanism of light emission due to exciton 
recombination in semiconductor CNTs has been proposed and experimentally verified^^i^i. 
A possibility of terahertz emission in CNTs imposed to transverse and axial electric field 
due electric-field induced heating of electron gas has been investigatecP^'^'^'^. Recently, 
the idea using kinetic energy of CNT-guided electron beam for stimulated emission of elec- 
tromagnetic waves in optical and terahertz ranges has been proposed^SiSSiSi. In the given 
paper we present a consistent theory of the effect. 

There is a wide family of devices utilizing interaction of electron beams with electromag- 
netic waves to produce electromagnetic radiation. Started by the invention of klystrons^^. 



this family embraces such well-known systems as traveling wave tubes (TWT) and backward 
wave oscillators (BWO)^, free electron lasers (FEL)^i^i2£i21^ etc. In systems of that kind, 
synchronous motion of electrons and electromagnetic wave modulates the electron beam and 
coherent radiation is produced by electron bunches. The radiation frequency is smoothly 
tunable due its dependence on the electron beam energy. Therefore, such type of emitters 
can operate in wide spectral range from microwave and infrared frequencies to vacuum ul- 
traviolet nowadays (e.g., VUV-FEL at DESY). Several projects aimed with the lasing in 
hard X-ray range have startecl2^>^. 

The synchronization of moving electrons and electromagnetic wave is attained ei- 
ther by slowing down the electromagnetic wave (Cherenkov, Smith-Purcell^ and quasi- 
Cherenkov^ radiation mechanisms) or by applying an external magnetic field, which is 
uniform in gyrotrons^i^ and spatially periodical in undulators^^. Besides, the oscillator- 
type mechanisnt^ is realized for electrons with discrete spectrum of transverse motion (for 
example, for electron channeling in crystals). The Cherenkov radiation is governed by the 
synchronization condition uj — ku = 0, where k is the wavevector and u is the charged 
particle (electron) velocity. In systems with external fields the synchronism condition is 
transformed to uj — ku — fi = with Vt as the electron oscillation frequency. In the oscillator 
regime VL is the transition frequency between electron levels^. 

For the coherent generation in the devices described, a high vacuum must be maintained in 
the region of the electron beam - electromagnetic wave interaction^. Otherwise, collisions 
of electrons with atoms move electrons out the synchronism and, consequently, lasing is 
not reached. From this point of view CNTs are unique objects since they exhibit ballistic 
electrical conduction at room temperature with mean free paths on the order of microns 
and even tens of microns^i^i^. Therefore, electrons can emit coherently from the whole 
CNT length which is typically 1-10 /im. Besides, single- and multi-walled carbon nanotube 
can carry a high current density of the order of 10^ — lO^'' A/cm^, see e.g. Refs.— i^"^. 
Lastly, metallic CNTs exhibit a strong, as large as 50-100 times, slowing down of surface 
electromagnetic waves^"^. Thus, a combination in CNTs of three key properties, 

(i) ballisticity of the electron flow over typical CNT length, 

(ii) extremely high current-carrying capacity, and 

(iii) strong slowing down of surface electromagnetic waves, 
allows proposing them as candidates for the development of nano-sized Chernekov-type 



emitters - nano-TWT, nano-BWO and nano-FEL. 

The remainder of the paper is organized as follows. In Sect. [TTl we derive dispersion 
equation for electromagnetic wave coupled with electron beam and discuss its solution in 
classical and quantum limits. A solution of the boundary- value problem for a finite-length 
CNT is presented in Sect. IIIII allowing evaluation of the absolute instability gain and the 
lasing threshold currents. Section IIVI contains physical analysis of the results obtained and 
numerical estimates for the gain and threshold currents. Concluding remarks are given in 
Sect. |Vl 

II. RADIATIVE INSTABILITY OF ELECTRON BEAM IN CARBON NAN- 
OTUBE 

A. Self consistent equation of motion for electromagnetic wave and electron beam 

Nanotubes - quasi-one-dimensional carbon macromolecules - are obtained by rolling up 
of graphene layer into a cylinder. The transformation can be performed in different man- 
ners classified by the dual index {ni,n2). The two integers rii and n2 represent the vector 
characterizing the way of turning, with rii = for zigzag CNTs, ni = n2 for armchair 
CNTs, and < rii 7^ ?t,2 for chiral CNTs. A nanotube can manifest either metallic or 
semiconductor properties, depending on its radius i?cn and the direction of rolling up. This 
correlation arises from the transverse quantization of charge carrier motion and is due to 
the quasi-one-dimensional topology of CNTs; for details see, e.g., Refs.-i^. 

Consider an electron beam moving in an isolated single-walled carbon nanotube oriented 
along the z-axis. The electron beam can be injected into the nanotube from outside by an 
external source or can be produced by applying voltage to some section of the nanotube. 
Accelerated by the voltage, electrons are injected into the working region. Independently 
on the origin of electrons, their motion in this region is assumed to be ballistic. 

As was mentioned in introduction, there is a certain analogy between a CNT guiding 
electron beam and macroscopic vacuum electron devices. The main (and obvious) distinction 
is the small cross-sectional radius of CNTs as compared to their macroscopic analogs. As a 
result, in CNTs spatial quantization of the electron motion comes into play and, therefore, 
classical models of the electron beam becomes inapplicable. The electron motion in CNTs 



is governed by quantum-mechanical equations. In this paper we shall consider the lasing 
effect when generated field is rather large, i.e., the condition 

E^Vh^(-Y (1) 



v c- 

is fulfilled^. In this case the electromagnetic wave has classical character and is described 
by the classical wave equation: 

VV-E(r,cu)-AE(r,a;) = ^^j(r,a;). (2) 

If the condition ([T]) does not hold, the number of photons per quantum level becomes too 
small to apply a classical approach and the electromagnetic field must be considered within 
the quantum electrodynamics. The quantum-electrodynamical consideration is of impor- 
tance on the initial stage of the instability development, when few photons participate in 
the process. We leave this stagefor further analysis focusing on the stage of highly devel- 
oped instability. Thus, in our model the electron motion is governed by the Schrodinger 
equation while the electromagnetic field is described by classical Maxwell equations. In the 
right-hand part of the field equation ([2]) the quantity j(r, uj) is the current density averaged 
over the quantum states of the electron beam. 

The current density in the working region is defined by the well-known equation^: 

j(r, t) = ^ {r{r, t)pV^(r, t) - (pV^*(r, t))V^(r, t)} 

-^|V^(r,t)|^A(r,t). (3) 

irieC 

Here p = —ihd/{dr) is the momentum operator and A(r,t) is the vector potential of elec- 
tromagnetic field. Further we neglect the Fermi law for the electron statistics. This is 
possible because the number of excited electrons per quantum level is found to be small 
even at superior current densities reachable in CNTs^'^'^. Indeed, number of levels in the 
interaction volume V is estimated as ~ Vp^/{2TTh)^, where p is a typical value of quasi- 
momentum of electrons in the beam. The number of electrons in this volume is ~ rieV, 
where He is the electron density. Then, the number of excited electrons per level is given 
by rje = {2TTh)^ne/p'^. At current density 10^ to 10^'^ A/cm^ and an excitation energy of the 
order of several electronvolts, we find rje ~ 10~^ — 10"'^. Therefore, the exchange interaction 
between electrons in the beam can be neglected. 



Let ip{v,t = 0) = ipn{'^) be the eigenfunction of an electron noninteracting with elec- 
tromagnetic wave and moving along the CNT. When the interaction is switched on the 
wavefunction is represented by the expansion 

ilj{r, t)=Y^ ai(t) exp{-ieit/h)iJi{r) (4) 

I 

over a complete set of the unperturbed eigenfunctions ipi{r) with ei as corresponding en- 
ergy eigenvalues. For further convenience, we rewrite the coefficients a/(t) as ai(t) = 
5in + 5ai (t), where 5in is the Kronecker symbol. Corrections ^aj" (t) are due to the electron- 
electromagnetic field interaction. Taking into account axial periodicity of the nanotube 
potential, the wavefunctions ipi{r) can be written in accordance with the Bloch theorem as 

iJi{r) = exp {ipiz/h} ^ bir exp {irz} M;^(r_L). (5) 

T 

Here pi is the axial projection of the quasi-momentum of l-th state, bir are constant coef- 
ficient, r = 2TTq/a are the reciprocal lattice constants, a is the CNT spatial period in the 
axial direction, uir {rj_) are functions dependent only on transverse coordinates, and q are 
integers. The term ^ 6;^ exp {irz} uir{rj_) is periodical in the z direction. 

T 

In linear approximation, the contribution to the electron current ([3]) originated from the 
electron-electromagnetic field interaction is described by the equation: 

5kir,t) = ^^{(5aS")*(t)exp[z(£, -£„)V^][V^r(r)p^n(r) - (M*(r))V^n(r)] 
^ I 

+ 6a^\t)exp[-z{ei-en)t/h][rn{r)pMr) - (pC(r))^Kr)]} 

^' |^„|'A(r,t). (6) 



mf,c 

Then, applying to Schrodinger equation standard perturbation-theory technique^ we obtain 
the equation describing the dynamics of the coefficients Sai{t): 

.^sr^ d6a,(t) , , , , . ,^. e 

X [A(r, t)p + pA(r, t)] ^„(r) exp(-i£„t/n) , (7) 

which is obtained by substitution of (j4]) into the Schrodinger equation and its subsequent 
linearization with respect to the electromagnetic field strength. The Fourier transform of 



Eq. ([7]) gives 

+pA(r,u; + ^^)|n). (8) 

Here we use the standard ket- and bra- notation of wavefunctions and matrix elements, 
10 — 4'i{^)- Only that terms are preserved in ([8]) which correspond to resonant interaction 
between electrons and electromagnetic field. Contribution of the last term in iQ is therefore 
neglected in ([8]). Performing the Fourier transform of Eq. ([6]) along the axial coordinate 
and time, we come to the k, w-space interaction-induced current density correction: 



bl,bnr [Ul, (Pn + t) + (p„ + t) U^,] Unr 
hu! + ei{pn- k) - EniPn) 

I Krhr' [Kr (Pn + t) + (p„ + t) M^J Uir' \ 

^ + £n{Pn) - £l (Pn + k) ) ' 

For convenience, we have introduced the vector form for the lattice constant r: t = re^, 
where e^ is the unit axial vector. The quasi-momentum operator in matrix elements is given 
by Pn = {'P±,Pn}, where axial components Pn are C-numbers and transverse components 
p_L are operators. These operators act only on the right- adjacent functions. Deriving (Q, we 
neglect the longitudinal component k of the electromagnetic wave vector in matrix elements 
since hk/pn <^ 1. Summation over the lattice constants r and r' is not independent: for 
every r in sum, the value of r' must be such that the values Pn + t — r' are in the first 
Brillouin zone. The coefficients Bni{k, rj_, uj) are given by 

Bniik,r±,uj) = ^bir'bl^ {unr\it>n + t) A{k,r±,uj) 

t't 
+A{k, r_L, Uj) {pn + t) I Ulr') . 

Substituting then ([9]) into (j2]) we come to a self-consistent field equation necessary for the 
further analysis. 

B. Dispersion equation for electromagnetic wave coupled with electron beam 



a 



Electromagnetic response properties of an isolated single-walled CNT was studied in Ref. 
on the base of a tight-binding microscopic model of the CNT conductivity and the effec- 



tive boundary conditions for electromagnetic field imposed on the CNT surface. A detailed 
analysis of the eigenwave problem has revealed propagation in CNTs strongly slowed down 
surface waves allowing the concept of nanotubes as surface-wave nanowaveguides. Consider- 
ing the electron beam as a perturbation, we can use the dispersion equation for the surface 
waves and the propagation constants obtained in Ref. |7| as a zero-order approximation. 
Then, the self-consistent field of the electromagnetic wave coupled with electron beam can 
be presented by the expansion 

A{k,r±,uj) = '^amik,uj)Amir±), (10) 

m 

where vector potentials Am{r±) correspond to the electromagnetic field eigenfunctions eval- 
uated in Ref.- and am{k,uj) are the coefficients to be found. Substitution of flTOj) . Qj and 
([6]) into ([2]) gives the system of equations for the electromagnetic field interacting with the 
electrons occupying n-th state: 
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bl,h^r [ul, (p« + r) + (p„ + t) u^^,] u^^ 

hw + ei{pn-k) -EniPn) 

I KAr' Kr (Pn + t) + {Pn + t) <,] Ui^' \ 

fkj + EniPn) - £1 {Pn + k) j' 

Here km are the wavenumbers corresponding to the physical system devoid electron beam. 
As one can see, deriving ( fTTl) we have proceeded from the single-electron dynamics to the 
dynamics of the electron beam: rig is the electron density. Multiplying left- and right-hand 
parts of Eq. flTTl) by A^(r_L) and utilizing the wavef unctions' orthogonality, we come to the 
dispersion equation as follows: 

, ,2 

^ r V"~^ , „(m) |2 



k-k ^ Vii?^^ 



^m"''e^ 



X 



+ 



-huj + SniPn) - £1 {Pn - k) fjhJ + SniPn) " ^i (p„ + k) 



(12) 



,(m) 



The upper index in Bl^^ relates the matrix element with the corresponding mode of the 



electromagnetic field A^ (r±), ^l = 2^J'Ke'^ne/me is the Langmuir frequency of the electron 
beam. 



The transcendent dispersion equation (IT2|) predicts the existence of a variety of branches 
of wavenumber k. Among them, the number of branches to be accounted for is defined by 
specific physical parameters of analyzed system. In the vicinity of a resonance, only terms 
corresponding to the resonant interaction, one or several (in the case of level degeneration), 
can be kept in the dispersion equation. If the difference between levels exceeds the linewidth, 
the only the resonant term is of importance. 

C. Classical and quantum limits in synchronism conditions 

Two terms in the right-hand part of Eq. flT2|) dictate two synchronism conditions corre- 
sponding to the resonant interaction between electron beam and electromagnetic wave: 

±hu! + EniPn) - El (,Pn ± k) = 0. (13) 

The signs " +" and " — " correspond to the absorption and the emission of photon by electron, 
respectively. Dependently on the relation between electron and photon energies, different 
interaction regimes are realized. As we restricted ourselves to the case when the photon 
momentum is much less than the electron one, the electron energy Ei {pn ± hk) can be 
presented by the truncated Taylor series as 

OEi(v ) 
Ei{pn ± hk) = Ei{pn) ± hk = Ei{pn) ±hkvi, 

OPn 



where vi is the electron group velocity. Then, denominators in (1121) can be represented by 

±hw + En{Pn) - ei{pn ± k) 

^±h{oo-kvi±n^i) + ^^k\ (14) 

z op^ 

The first term in the right-hand part of this equation is analogous to the standard term 
uo — ku ±VL in the synchronism condition'^^ . The only difference is that the velocity of 
free electrons is replaced by the group velocity of quasi-electrons vi and the undulation 
frequency is replaced by the transition frequency fi„; = [Enipn) — Ei{j)n)]/h between CNT 
energy bands. The last term in flT4l) originates from the quantum recoil of electron during 
emission (absorption) of photon and induces a red (blue) shift of the transition frequency. 
This term is inversely proportional to the electron effective mass (second derivative of the 



energy). Let / = s be an electron level corresponding to the resonant interaction. Then, 
within the approximation stated, the dispersion equation takes the form as follows: 



ijj — kv^ — ———-^ — Q 



ns 



2 dpi, 
where 

~ 8m A;' c^ ' "■^ ' ' ™ ~ rtei^AimJ- 
In the case of intraband transitions Qns = and Eq. ( TT5l) takes the form of the dispersion 
equation for the instability with the recoil accounted for—. 

Depend on ratio between the radiation linewidth and the recoil-induced detuning, two 
different generation regimes are realized. In the low-gain limit^ the spontaneous emission 
linewidth can be estimated as Auj/uj ~ c/{ujL), where L is the interaction length. If the 
linewidth exceeds the recoil energy, the recoil term in the denominator of ( IT5l) can be ne- 
glected and the classical interaction regime is realized. The dispersion equation in that case 
takes the traditional form of the second-order Cherenkov resonance: 

dpi {uj - kv,Y 

The spatial increment of the instability k" = Im(A;) can be estimated using the method of 
weakly coupled modes^. According to this method, interaction between the electromagnetic 
wave and the electron beam is essential only in the vicinity of the point {uq, ko = ujq/vs) 
where the dispersion curves of noninteracting modes, u — kvg = and k{ij) = km{uj), are 
crossed. Then km is represented by the expansion 

dkmiuj) 



km{uj) = ko + 



du 



{oo - loo) . (17) 



Substitution of this expansion and k = ko + Ak into ( fT6l) results in a third-order alge- 
braic equation with respect to Ak. From this equation, the instability spatial increment is 
estimated at the frequency u = ujq as 



s/3 



lAri ^ 



,(m) O En K, 



1/3 



) 



(18) 



where Ak" = Im.{Ak). Since 6„„ ~ ng, the increment is found to be the 3-rd root of the 
electron density. Such a dependence is typical for the Compton-type radiative instability^^. 
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In the opposite case, when the hnewidth is less then the difference between the emission 
and the absorption frequencies, we fall into regime of the strong quantum recoil impact. In 
this case, only the term corresponding to the emission survives in the dispersion equation 
flT^ . which therefore is reduced to 



k-k^ = ^ ^^_ . (19) 

As a result, the instability increment is given by 



^ dPn 



\Ak" 



Am) 



hv„ 



1/2 

(20) 



i.e., turns out to be proportional to the square root of the electron density. 

Below we present a detail discussion of the different generation regimes and give some 
numerical estimates of physical parameters corresponding to these regimes. 

III. STARTING CURRENTS AND THEIR DEPENDENCE ON THE NAN- 
OTUBE LENGTH 

A. Boundary conditions for a finite length nanotube 

In sections IIIBI and IllCt dispersion equations have been derived providing us with 
wavenumber eigenvalues in an infinite-length CNT guiding electron beam. As a next step, 
the system must be imposed by edge conditions accounting for the finite length of the inter- 
action zone. These conditions are stated as the requirement to perturbations of the electron 
and current densities, generated by the electron beam - electromagnetic wave interaction, 
to be zero at the input of the working zone, i.e. 

5n,{z = Q) = 53n{z = Q) = Q . (21) 

The condition that the tangential electric field component and the axial component of the 
magnetic field be continuous on the CNT surface yields additional boundary condition. We 
write it in the simplified form^- as 

E{z = 0) = aE{z = L) , (22) 



11 



where a is the reflection coeflicient of electromagnetic fleld from the working zone boundaries. 
The field distribution in a finite-length system consisting of several parts can be found by 
solving electrodynamical problem in each region separately and then joining the solutions by 
means of boundary conditions. In the interaction region, the electromagnetic field is given 
by 

N 

E(z) ~^Ciexp(zfc(^)z), (23) 

i=\ 

where the summation is performed over all electromagnetic modes in CNT; the wavenum- 
bers fc^*'' are determined by corresponding dispersion equations. Note that the refiection of 
electromagnetic waves from boundaries back into the working zone creates positive feedback 
in the system and thus allows accumulation of the electromagnetic energy and provides an 
oscillator regime. 

B. Starting current at a large quantum recoil 

In the quantum interaction regime, when the quantum recoil exceeds the linewidth, the 
instability is described by the quadratic dispersion equation flT^ with solutions k^^^ and hS'^^ . 
Consequently, the electric field and the perturbation of the current density in the working 
zone are given by 

E ~ ciexp(ifc(^)z) + C2exp(zA;(^)2) , (24) 

5j„ ^ ^exp(^A;(i);2) + ^ exp(^A;(2)^) . (25) 

The coefficients 

K2 = 1 - -k^'''^ + ^^ k^''"^" (26) 

introduce deviations of the wavenumbers k'^^^ and /c*^^-* from the synchronism, and the co- 
efficients Cj are determined from the boundary conditions as was discussed in Sect. IIII A[ 
Using the boundary condition fl2T|) and (!22|) . we arrive at the linear system for Cj as follows: 

ci + C2 = a[ci exp(zA;*^^)L) + C2exp(iA;(^)L)] , 

ci C2 (27) 

— + — =0. 

Nontrivial solution of this system is determined by the equation 

(5i[l-aexp(zA;WL)] - 52[l - ae^^{ik^'^^ L)] =0. (28) 

12 



A current density satisfying equation (128!) is the threshold current density of the generation. 
To evaluate this quantity, characteristic equation fl28|) must be solved together with equation 
(TT^ . Substituting the roots 

fc(^'^) = km,. + , '"" ,.. (29) 

of the dispersion equation ( TT9l) . with /cch extracted from the synchronism condition uj — 
kchVn + (^^ch/2)<9^^n/c^Pn = 0' i'^to (l28l) and solving the resulting equation with respect to 
the current density, we obtain 

iXm) ■ 2 

-^U—^ = l-\a\+Lkl, (30) 

where 

^ = r " ^™-n + ^^J - (31) 

is the dimensionless off-synchronism parameter. 

Physically, Eq. (130|) establishes the energy balance in the working zone. Its left-hand part 
determines the radiation production which is therefore proportional to the electron density 
rie and to the squared interaction length. The factor sin^ x/x^ determines the so called gain 
curve — the gain dependence on the off-synchronism parameter x. In the case considered the 
gain curve is symmetrical with respect to x = and is maximal at zero deviation x. Further 
we compare this result with the classical case of small recoil and demonstrate significant 
difference in the behavior of gain curves. The term 1 — a in the right-hand part of fl30|) 
corresponds to the radiation leakage through the boundaries of the interaction zone while 
the last term specifies the radiation absorption by nanotube. 

The energy balance equation (l30ll allows the evaluation of the threshold current density. If 
the current density in CNT exceeds the threshold value, the generation process is developed. 
The characteristic time of the instability development is inversely proportional to the absolute 
instability increment uj" = Im(ci;), which is derived by solving the generation equation (!28|) 
with respect to uj{k). In the lowgain regime^, which implies the conditions \Ak"\L <^ 1 and 
1 — a ^ 1 to be fulfilled, the increment is given by: 

dkr, 



< 



did 



'™L^-i^-A;M. (32) 



In the linear stage of the radiative instability development, the electromagnetic field grows 
with time as exp(ix;^t). 

13 



C. Starting current in the classical regime of interaction 

In the case when quantum recoil can be neglected, the dispersion equation (fT6l) gives 
three roots 



dpi {UJ - Vnk'j' 

k^'''^ = A:ch ± -Jbt^^jr^. (33) 

and, consequently, electromagnetic field in the interaction region is given by Eq. (1231) with 
A^ = 3. Correspondingly, perturbations of the electron and the current densities in the beam 
are written as 

3 3 



j=l « j=l 



where deviations z/j are given by fl26|) with the last term omitted, i.e., z/j = 1 — k^'^^Vn/uj. 
Then, by analogy with the previous section, we obtain the linear system 



Ci + C2 + C3 = a [ci exp(zA;(^)L) 

+C2exp(zfc(^)L) + C3exp(zfc(^)L)] 



Ui U2 U-i 

ui vi vi 

and corresponding generation equation 



(35) 



-? 



Z/2-Z/3) [l-aexp(2fc«L)] 
-vl{y^ - z/a) [1 - aexp(zfc(2)L)] 
+z/|(z/i - z/2) [1 - aexp(iA;(3)L)] =0. (36) 

This equation we solve in the low-gain limit, which is determined by the condition k'^L < 1. 
The curve depicted in Fig. [T] divides out areas of parameters corresponding to low- and 
high-gain regimes, respectively. Then, solutions of Eq. fl36l) - the threshold current and the 
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low-gain region 



CNT length, /an 



FIG. 1: Demarcation between low- and high- gain regimes of the generation equation (j36p . 



temporal instability increment - are given by 



dm) 



d'^e. 



vl dpi 



kU 



X cos X — sm X 



UJ„ 



Okrri 



X'^ 

dm) 
Jnn 



1- \a\ +Lk'' 



(37) 



vl dp. 



d^En j2 ^ "^os X — sin x 

-Lj 



X"" 



\a\ 



L 



+ K 



(38) 



with the parameter x defined by Eq. fl3T]) . As follows from the balance equations fl30|) 
and (1371) . in the quantum interaction regime the radiation production per unit length is 
characterized by the linear dependence on L, while this dependence becomes quadratic 
in the classical regime. Besides, the gain curves display distinctive behavior in these two 
cases. As different from the quantum interaction regime, in the classical limit the gain curve 
has asymmetrical character^^ due to the interference of absorption and emission processes 
separated in this case by a frequency gap narrower then the linewidth. As a result, sign 
of the absolute instability increment depends on the sign of the synchronism detuning. At 
positive detuning the system is closer to the absorption frequency while negative detuning 
moves the system to the emission frequency. 

Qualitatively, the classical and quantum interaction regimes are divided by the demarca- 
tion line depicted in Fig. [2l In the area above the line the quantum recoil at the generation 
must be taken into account while in the area below the line this effect can be ignored. The 
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large recoil impact 



small recoil impact 



CNT length, //m 

FIG. 2: The curve divides the regions of parameters with small and large impact of the quantum 
recoil on the generation. The curve has been obtained for the low-gain regime. 

line course can easily be explained just by the increase of the photon energy with frequency. 
In addition, the increase of the generation length L leads to narrowing of the gain line and, 
as a result, the quantum recoil conies into play at smaller frequencies. 



D. The role of electron spread 

If electrons in the beam are distributed over a large number of energy levels and energy 
spread significantly exceeds the gap between emission and absorption lines, the total current 
is obtained by summation over this distribution. The generalization of Eq. flT6|) on this case 
is obvious: 



k - k^co) = -fei'::) / dvfiv) 



d'Sr. 



e 



(39) 



dpi {io -vkf 
This equation can be rewritten in the form, conventional in plasma physics^: 

Deriving (!39l) we assumed the dependence of the normalized distribution function 
(/ f{'^)dv = 1) on the group velocity to be narrower than corresponding dependences of 
the matrix element form and second derivative of the energy d'^en/dpl. Then, considering 
the group velocity spread exceeding the spontaneous emission linewidth, uAv/c ^ c/^uL), 
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in (139|) we can make use the standard representation 

= V - m5 (to - vk) . (40) 

The principal value of the integral determines the real- valued component which is out of our 
interest . 

If resonant interaction between electron beam and electromagnetic field occurs in the 
region of the negative derivative of the distribution function, i.e. df{v)/dv < 0, then k" > 
and the generation process is not developed (we choose the exp{ikL) dependence). This 
is because the majority of electrons in that case have velocities smaller then the resonant 
velocity and therefore they absorb the electromagnetic wave energy. Such a situation takes 
place in equilibrium, when the number of particles occupying energy level grows less with 
the level energy increase. In such a system, an initial perturbation attenuates. This process 
is commonly known as the Landau attenuation. 

If the resonance is in the region with positive derivative df{v)/dv > 0, the radiative 
instability is possible and obeys the condition 






d'sn df{v) 



<0, (41) 

v=uj/k 



"^ "" dpi dv 

which originates from the requirement k" < and from Eqs. fl39|) and fHOj) . The condition 
(14T|) expresses the excess of emission over absorption. As one can see, the emission per unit 
length does not depend on the interaction length. 

The imaginary part of the wavenumber k describes the asymptotic exponential behavior 
of the electromagnetic field in a continuous medium. To reach generation in a finite region, a 
corresponding boundary conditions must be imposed. At a large spread, when the resonant 
term in (l39l) can be presented by (l40l) . the dispersion equation (139!) has the only root. Using 
fl22|) we arrive at the relation ci = aexp [ik^^^LjCi , which dictates the generation equation 
as 

l-aexp(ifcWL) =0. (42) 

For the Cherenkov radiation mechanism, solution of (H2]) leads to the equations as follows 



17 



for the threshold current density and the absolute instability increment: 

lal + Lk'' , 






dpi 



dv 



(jJrr 



dk 



du 

1 — la 



TT. 



= 1 

v=uj/k 

d'Sn df{v) 



(43) 



5M 



Qp2 Qy 



v=uj/k 



L 



k'L 



(44) 



Equation ( l44l) shows that the production of stimulated radiation in the case of large spread 
is defined by the spread and falls down with its increase. The line dividing the range 
of parameters into two domains, with weak and strong influence of the energy spread, is 
depicted in Fig. [31 With the CNT length increase the role of the spread also rises due to 
the gain line narrowing. 



large spread influence 



4 6 

CNT Icngtth, /m 



FIG. 3: The regions of parameters with small (below curve) and strong (above the curve) influence 
of the electron beam spread. 



The extension of the obtained generation conditions to the case of interband transitions 
(i.e., to the undulator regime) is obvious and, in accordance with the dispersion equation 
flTSj) for the undulator regime, is achieved by the substitution uj — kvn —>■ uj — kvg — ^ns in the 
off-synchronism parameter (!3T|) . and the substitution kk'^d'^En/dp'^ -^ kk'^d'^Eg/dp'^ — 2Qns 
in expressions for the threshold current fl5Ul) . (157]) and fH51) . and for the absolute instability 
increment (BSD, (1551) and 
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The analytics presented in this section imphes fulfillment of several simplifying approx- 
imations: smallness of the photon momentum hk/p <^ 1, small or large infiuence of the 
electron recoil on the emission (absorption), small or large electron spread). Obviously, the 
analytical approaches do not work in intermediate cases; equation flT^ supplemented by cor- 
responding boundary conditions requires numerical integration. The number of roots of the 
dispersion equation to be accounted for and corresponding number of boundary conditions 
to be imposed is dictated by concrete physical parameters of the system being considered. 

IV. PHYSICAL ANALYSIS AND NUMERICAL ESTIMATES 

In Sect. IIIAl it has been stated that classical treatment of electromagnetic field is valid if 
the field strength Ee amounts to a certain sufficiently large value. This value is determined 
by the condition imposed on the number of photons per energy level to exceed unity^. At 
the initial stage of the instability development, with less then one photon per energy level, 
the photon dynamics is described within the quantum electrodynamics formalism. 

Usually, the number of photons per energy level is given by riph {c/uj) , where Uph is the 
photon number per unit volume while the quantity {uj/cY determines the number of photon 
levels lying below the energy hw. As different from that, in the case of high-coherent laser 
radiation the radiation is concentrated in a narrow spectral range Acj ~ c/ L. As a result, 
the parameter defining the possibility of classical consideration of electromagnetic waves - 
the number of photons per energy level - is derived as the density of the beam's kinetic 
energy converted to electromagnetic field divided by the photon energy and the number of 
levels below t\w. The ratio is found to be 

i /c\2 /c^cmc^(7- 1) 
~ r]ph— [-] L- 



ve \uJ u hjjj 

where j is the current density, v is the electron velocity, rjph is the efficiency of the transfer 
of electron kinetic energy to electromagnetic field. For infrared photons and electrons of 
several electron- Volt energy and ~ 10 /xm length nanotube, the photon number per energy 
level exceeds unity (i.e., the classical treatment is possible) if rjph > 10^^. Since the initial 
stage of the instability development is beyond the scope of our paper, the parameter rjph can 
be estimated from the relation rjph ~ l/{kL) ~ 0.02, which corresponds to so-called nonlin- 
ear saturation regime^ and determines the electron beam energy conversion in saturation. 
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Therefore, generation threshold and nonhnear stage of the instabihty development can be 
considered classically. 

A simplest way to realize nanoFEL in carbon nanotube is to inject into it a high energy 



external electron beam. Since the velocity of free electron is f (cm/s) = 5.7 x 10"^ ^^/ip{eV) , in 
order to accelerate electrons up to velocities providing the synchronism regime (with 50-100 
times wave slowing down predicted in Ref. uj), it is necessary to apply voltage of (y9 ~ 7 eV. 
If the CNT diameter is such that its product with the electron transversal momentum is 
p±D/h ~ 10 — 100, the electron motion can be treated as classical. In that case, the term 
in the right part of the dispersion equation ( TT6l) can be modified in the following way 



j,(m) 



d^Er 



k' 



ve 



~ UJ, 



(45) 



dpi {u-vnkf ' "-^k'^c^^uj-vkf 

where v is the classical electron velocity and e is polarization vector for the electromagnetic 
mode considered. This simplification, after substitution of f H5]) into Eqs. (1571) and (I5S1) . 
allows us estimate the threshold current required to start the generation process and the 
instability increment, respectively. The dependences of these quantities on the CNT length 
are depicted in Figs. H] and [51 Calculations have been done for 1 /im radiation wavelength 
and for the reflection coefficient from the working zone boundaries a = 0.99. Generation in 
the terahertz range would require higher current density. It follows from Fig. [Hthat the gain 
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FIG. 4: The dependence of threshold current density on nanotube length. 



for CNT is extremely large as compared with macroscopic electronic devices. For chosen 



20 



a loo- 



se 



^ 











a: j= lOxloVcm^ 
b: j= 7 X loVcm^ 
c: j= 5xl0'A/cm^ 


a / 


■ 


^ 


y^ b 


_ ^ 


^ 


^ 15 
CNT lenght, /an 



FIG. 5: Instability increment vs nanotube length at different electron current densities. 



parameters, the generation development starts when the CNT length is about 6 microns 
or larger, what is technologically routine range. Therefore, our calculations demonstrate 
that the development of CNT-based nanoFEL is already possible at the current stage of 
nanotechnology. The characteristic time of the instability evolution is inversely proportional 
to the instability increment and for 10 /im nanotube is a fraction of nanosecond. 

A positive feedback is required for the realization of oscillator regime; reflection from 
CNT ends^- can serve as a possible mechanism of the feedback. The reflection can by 
intensified by variation of the CNT generic parameters, proper selection of surrounding 
medium and using other methods commonly applied in laser physics and electronics. An 
alternative mechanism providing the feedback is excitation of backward modes propagating 
oppositely to electron flow. The backward modes are possible because CNTs are periodic 
along their axis and, consequently, their eigenmodes are Bloch modes containing waves with 
both positive and negative phase velocities. As a result, there exist Bloch modes with group 
velocity directed oppositely to the electron velocity - the backward modes^l. One of the 
waves of the backward mode having a positive phase velocity can be synchronized with the 
electron flow. In this case the positive feedback is provided automatically. 

The instability process is developed only if the electron free-path length is comparable or 
even exceeds the working zone length, i.e., the electron motion is ballistic within the zone. 
Otherwise, random collisions of electrons cause a phase shift which prevents the electron 
flow bunching and brakes the radiation coherence. As was mentioned above, in metallic 



single-walled CNTs the free-path length is about several microns 



47.48.49 



A longer ballisticity 
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area can be provided by proper external conditions. For example, in a regular array of 
oriented nanotubes the suppression of electron collisions with atoms can be achieved using 
the properties of the electron diffraction in periodical structures. In a densely packed array 
of CNTs — CNT bundle — nanotubes form a lattice with the distance between CNTs' 
axes 2i?cn + d, where d ^ 3.2 A is the interlayer distance in graphite. Correspondingly, 
the reciprocal lattice vector in such a lattice has the value h = 2'7r/(2i?cn + d). From the 
principle of uncertainty we can estimate the transverse component of the momentum by 
p_L/h ~ 2iT/2Rcn- Obviously, the Bragg condition |p^ + h\ ^ |p_|_| can be fulfilled for a large 
portion of electrons passing the bundle and six-wave diffraction^^ can be realized. Owing to 
the diffraction, electrons are concentrated in domains free of atoms and, therefore, scattering 
is weak for such electrons. Analogous situation meet in the Bormann effect^ for hard X-rays 
passing through a crystal. Owing to this effect, a significant increase of the photon free path 
is observed. 

Even if generation conditions are provided by the use of external electron beam, the 
idea to exploit intrinsic electrons of CNTs looks very attractive because it would solve the 
dramatic problem to focus an external electron beam into a spot of the CNT diameter size. 
Typical velocity of 7r-electrons excited to energy of several electronvolts is about'^ 10^ cm/s. 
For such electrons, the synchronism condition requires the electromagnetic wave slowing 
down as large as 300 times, which is much larger than the theoretical estimate^ gives for 
CNTs. 

In such a situation, special configurations providing higher group velocity are extremely 
desirable; otherwise, stronger excitation of electrons is necessary to fulfil the Cherenkov 
generation condition. Fortunately, as compared to vacuum electronic devices, stimulated 
emission in CNTs features a set of new promising properties. In macroscopic Cherenkov 
FELs the electron energy ordinarily rises with the electron velocity and, in nonrelativistic 
regime, quadratically depends on the momentum (and velocity). As a result, the only way to 
reach the synchronism condition in that case is to increase the electron beam energy. For a 
collective (quasi-) electrons in CNT such is not the case. Indeed, the electron group velocity, 
which is analog of the velocity for quasi-particles, is determined by the properties of the whole 
system and may demonstrate nontrivial dependence on the quasi-momentum. Locally, the 
quasi-particle velocity may recede as energy rises. Correspondingly, local maxima of the 
group velocity may appear. If one seeks the synchronism condition for a low-energy quasi- 
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particle, it is advantageous to choose parameters in the vicinity of the group velocity local 
maxima. It allows attaining the synchronism in a relatively low accelerating potential and, 
therefore, significantly reduces the CNT energy load. 

Let us exemplify the statement considering an isolated straight (g, q) armchair CNT. The 
dispersion law of 7r-electrons in such a CNT is given by^ 



£i{p) = =F7o [1 ± 4cos(7r//g) cos(ap) 
+4cos^(ap)]^''\ 



(46) 



where 70 — 2.7 eV is the overlap integral, / 



2g, a = V3b/2h, 6 = 1.42 A is the 



interatomic distance in graphite. The upper and lower signs refer to the conduction and 
valence bands, respectively. The group velocity corresponding to this law is 



Vi = =f27oO' sin (pa) 

=F cos (vr//g) — 2 cos (pa) 
[1 ± 4 cos (vrZ/g) cos (pa) + 4 cos^ (po-)] 

Calculations of the energy and the group velocity for (10,10) nanotube by Eqs. 



(47) 
and dm) 



^,(PJ' eV a) 




u,(p) X 10 , cm/s b) 




FIG. 6: Energy (a) and group velocity (b) vs quasi-momentuni for (10,10) armchair nanotube. 
Numbers near curves are mode numbers I. 

are presented in Fig. El The curves in two figures can easily be correlated: The larger slope 
of the dispersion curve the larger the group velocity. Typical velocity of vr-electrons excited 
to energy of several electronvolts is about^ 10^ cm/s. For such electrons, the synchronism 
condition requires the electromagnetic wave slowing down as large as 300 times. A proper 
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choice of the excited state in the vicinity of the group velocity local niaxiniums allows 
essential weakening this restriction. 

The region in the vicinity of the group velocity extrenium is also attractive because of the 
weak velocity dependence on the quasi- momentum. As a result, in this region irradiation 
of photon gets the electron only slightly out the synchronism condition keeping high the 
probability to emit next photon. Due to that, the radiation effectiveness grows in the 
vicinity of the group velocity extremum. An additional advantage of the local maximum 
in the group velocity is the smaller negative influence of the beam energy spread on the 
generation effectiveness. Indeed, in the vicinity of the group velocity extremum the Taylor 
expansion of the energy does not contain linear quasi-momentum terms. As a result, a larger 
number of particles in a spreaded beam appears to be synchronized with electromagnetic 
wave. This effect is characteristic for quasi-particles and fully absent for free electrons. 

The effect of radiation instability in nanotube can be controlled by the variation of the 
electron effective mass. The smaller the mass the more responsive the electron is to pertur- 
bation, and the more likely an electron beam bunching. This means a faster development 
of the instability. The reciprocal electron effective mass is given by the quantity (9^£:„/9p^; 
therefore, the increase of the instability increment as the effective mass grows smaller follows 
immediately from the dispersion equations (fT6l) and (139!) . which involve the reciprocal mass. 

One more mechanism, which does not require large wave slowing down, is exploiting 
electron interband transitions. In this case, as follows from ( ITSl) . the resonance condition is 
ijj — Vgk = Qns (we suppose that transition frequency exceeds the term related to the recoil) 
and the radiation frequency can vary from infrared to ultraviolet. For interband transitions, 
single-particle spontaneous emission of electron (positron) beams emerging from outside 
into nanotube was considered by Artru et al.—. 

To weaken the requirement imposed on the electromagnetic wave to be slowed down to 
the electron velocity one can utilize the photon diffraction on a periodic lattice of carbon 
atoms in a nanotube. Resonance interaction takes place for harmonics corresponding to the 
reciprocal vector r satisfying the condition u — Vn{k + r) = 0. Then, taking into account 
the condition f„/c <^ 1, one can obtain: 

^r = . ^. 48 

Here riref = kc/u is the effective refractive index of corresponding mode. The spatial period 
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of a nanotube varies in wide range. For zigzag and armchair nanotubes it equals to 2.49 
A, while for chiral nanotubes the translation period achieves 10 nm and more depending 
on the nanotube indices. As a result, the generated wavelength varies from ultraviolet 
(for armchair and zigzag CNTs) to infrared range for nanotubes with translation period 
~ 2v37ri?cn- 

V. CONCLUSION 

In the present paper, aiming at the development of the physical basis of a new class of 
nano-sized light sources, we have investigated theoretically a recently proposed mechanism 
of the generation of stimulated electromagnetic radiation by electron beam in carbon nan- 
otubes. The basic idea exploits an analogy between CNTs and macroscopic electron devices 
and utilizes the effect of wave slowing down in waveguides. Three basic properties of carbon 
nanotubes, the strong slowing down of surface electromagnetic waves, the ballisticity of the 
electron motion over typical CNT length, and the extremely high electron current density 
reachable in CNTs, allow proposing them as candidates for the development of nano-scale 
Chernekov-type emitters for a wide frequency range from terahertz to optical. The thresh- 
old conditions evaluated from the theoretical model demonstrate that the development of 
CNT-based nanoFEL is already feasible at realistic present-day parameters of nanotubes. 
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